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Neutrons stars lighter than the Sun are basically composed of nuclear matter of density up to 
around twice normal nuclear density. In our recent analyses, we showed that possible simultaneous 
observations of masses and radii of such neutron stars could constrain ri = (KoL 2 ) 1 ^ 3 , a combination 
of the incompressibility of symmetric nuclear matter Ko and the density derivative of the nuclear 
symmetry energy L that characterizes the theoretical mass-radius relation. In this paper, we focus 
on the mass-radius constraint of the X-ray burster 4U 1724-307 given by Suleimanov et al. ,lj. We 
therefrom obtain the constraint that ri should be larger than around 130 MeV, which in turn leads 
to L larger than around 110, 98, 89, and 78 MeV for Ko = 180, 230, 280, and 360 MeV. Such a 
constraint on L is more or less consistent with that obtained from the frequencies of quasi-periodic 
oscillations in giant flares observed in soft-gamma repeaters. 

PACS numbers: 04.40.Dg, 21.65.Ef 


Introduction. Neutron stars, stellar remnants of su¬ 
pernova explosions at the end of massive stars, are con¬ 
sidered to be composed of matter in extreme condi¬ 
tions, namely, ultra-high density and large neutron ex¬ 
cess. Since the temperature of the matter is generally 
very low compared with the typical neutron Fermi tem¬ 
perature, it is extremely difficult to examine the equilib¬ 
rium properties of such dense cold matter in the labora¬ 
tory, although highly energetic heavy-ion collisions could 
create hot dense matter as encountered in protoneutron 
stars. Theoretically, on the other hand, the equation of 
state (EOS) for matter in neutron stars, hereafter re¬ 
ferred to as neutron star matter, remains to be deter¬ 
mined, particularly above normal nuclear density, p o. In¬ 
versely, neutron stars could be a suitable laboratory to 
probe the properties of cold dense matter. For example, 
observations of masses and radii of neutron stars would 
help us to constrain the EOS of neutron star matter. In 
fact, recent discoveries of neutron stars with about two 
solar mass ( Mp ) play a role in ruling out various soft 
EOS models 0, y|- Furthermore, estimates of radiation 
radii of neutron stars have been made via observations 
of thermonuclear X-ray bursts and thermal spectra from 
low-mass X-ray binaries , which could also give us a 
significant constraint on the EOS. Additionally, oscilla¬ 
tion spectra radiated from a specific kind of neutron stars 
are another observable information to see stellar proper¬ 
ties, such as masses, radii, the EOS, rotations, and mag¬ 
netic fields (e.g., [8t-ll5| ). This unique technique is known 
as neutron star asteroseismology. Although observational 
evidences for neutron star oscillations are extremely lim- 
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ited, quasi-periodic oscillations discovered in the after¬ 
glow of giant flare phenomena observed from soft-gamma 
repeaters [l 6 [ are considered to be strongly associated 
with oscillations of whatever portion of neutron stars. 
Through these observations, possible constraints on the 
stellar properties, particularly in the crustal region, are 
discussed jl7l - [23j . 

Since the details of neutron star structure obviously 
depend on the still uncertain EOS of neutron star mat¬ 
ter, they have yet to be clarified. It is generally consid¬ 
ered [ 2 J] that, under a liquid metallic ocean close to the 
surface, neutron-rich nuclei form a Coulomb lattice in a 
sea of electrons and, if any, dripped neutrons. Because 
of the crystalline structure, the corresponding region is 
called a crust. As the matter density increases up to a 
value close to po, it is considered that such nuclei begin 
to melt into uniform matter, which consists mainly of 
a core of the neutron star. Furthermore, non-nucleonic 
components such as hyperons and quarks might appear 
for a still higher density region inside the core, depending 
on the model for neutron star matter [25j |. In addition 
to the possibility that such non-nucleonic components 
appear, it is also suggested that the uncertainty from 
three-neutron interactions in the EOS for pure neutron 
matter comes into play for the same region [Hj. On the 
other hand, neutron star matter of density below about 
2po is relatively easier to be constrained from terrestrial 
nuclear experiments. This is why we will focus particu¬ 
larly on low-mass neutron stars that have central density 
p c lower than 2po. We remark that we succeeded in con¬ 
structing theoretical mass and radius formulae for such 
low-mass neutron stars, which are written as a function 
of p c and 77 , a combination of the EOS parameters that 
characterize the nuclear saturation properties [27] . 

In this paper, we systematically examine the p de- 
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pendence of the mass-radius relation of low-mass neu¬ 
tron stars using more than 200 phenomenological EOS 
models (HI that are constructed in such a way as to 
reproduce empirical masses and radii of stable nuclei. 
Then, by comparing the obtained mass-radius relation 
with available neutron star observations, we give possible 
constraints on r). In particular, for this purpose, we fo¬ 
cus on constraints on the mass-radius relation of neutron 
stars that were derived by Suleimanov et al. [l| from the 
observed cooling phases of the X-ray burster 4U 1724-307 
located in the globular cluster Terzan 2 via different at¬ 
mosphere models. This is because unlike other studies to 
make a constraint on the mass-radius relation via the ob¬ 
served thermal emission from neutron stars, Suleimanov 
et al. obtained such constraints by using the whole cool¬ 
ing track and checking the consistency with the theoret¬ 
ical prediction of neutron star cooling evolution, which 
enables us to minimize the theoretical uncertainties in 
the atmosphere model during the burst phenomena j 2 t| . 
As will be shown below, the resultant constraint on the 
density dependence of the symmetry energy is consistent 
with the known constraints from the quasi-periodic oscil¬ 
lations observed in giant flares of soft gamma repeaters 


EOS parameters. We begin with an expression for the 
EOS of uniform nuclear matter at zero temperature. The 
bulk energy per nucleon, w, of this matter can be gener¬ 
ally expanded around the saturation point of symmetric 
nuclear matter as a function of the nucleon number den¬ 
sity, n b, and neutron excess, a, as jJlO 


w = wo + -^2 (m>-n 0 ) 2 + 

18 u.q 


So + o—( n b - n 0 ) 
3n 0 


a 2 , ( 1 ) 


where a is defined as a = (n n — n p )/no with the neu¬ 
tron and proton number densities, n n and n p . That is, 
the case of a = 0 corresponds to symmetric nuclear mat¬ 
ter, while the case of a = 1 corresponds to pure neutron 
matter. The parameters wq, no, and Ko, which charac¬ 
terize this expansion, denote the saturation energy, sat¬ 
uration density, and incompressibility of symmetric nu¬ 
clear matter, respectively. On the other hand, So and L 
are the parameters associated with the symmetry energy 
coefficient, i.e., So is the symmetry energy coefficient at 
rib = no, and L is the density dependence of the sym¬ 
metry energy around rib = n-o- Note that among the five 
parameters in Eq. ©, wq, no, and So can be relatively 
easier to constrain from empirical masses and radii of sta¬ 
ble nuclei, while the remaining two parameters, Kq and 
L, are more difficult to fix j3l|. Thus, we particularly fo¬ 
cus this paper on the parameters Kq and L. We remark 
that many EOSs of nuclear matter have been proposed 
so far, which have various values of Kq and L, while hav¬ 
ing reasonable values of wo, no, and So (e.g., M). 
Although it may well be difficult to precisely describe 
the mass-radius relation of neutron stars by taking Kq 
and L alone as free parameters, these two parameters 
are expected to mainly control the stiffness of the EOS of 
neutron-rich nuclear matter near po and hence the struc¬ 



FIG. 1: (Color online) y as a function of L. The dots are 
taken from the 247 OI-EOSs, while the solid line denotes the 
fitting in a functional form of L 2 ^ 3 . 


ture of at least an outer part of neutron stars. In fact, we 
succeeded in finding a suitable combination of Ko and L, 
namely, 77 = (KqL 2 ) 1 / 3 , that well characterizes the struc¬ 
ture of low-mass neutron stars 27] in the sense that the 
mass-radius relation changes smoothly with 77 (see Fig. 
0 ). 

Now, in order to cover a wide range of 77, we consider 
the phenomenological EOSs of neutron star matter based 
on the simplified version of the Thomas-Fermi method 
that allows for the bulk, gradient, and Coulomb energies 
El- These EOS models were systematically obtained 
from the energy of uniform nuclear matter, which, in the 
limit of rib — >■ no and a —> 0, reduces to Eq. m with var¬ 
ious values of y = -~K 0 So/{3noL) and K 0 . In fact, the 
most relevant values of Wq, no, and So were determined 
together with that of the gradient energy coefficient for 
given y and Ko by fitting masses and charge radii of sta¬ 
ble nuclei obtained from the optimal nucleon distribution 
to the empirical ones [33]. We remark that y corresponds 
to the gradient of the saturation line near a = 0. Finally, 
the crustal EOS was obtained for various sets of ( L, K 0 ) 
[32 ] by extending the Thomas-Fermi method to several 
shapes of nuclei in a lattice within a Wigner-Seitz approx¬ 
imation J33j]. Hereafter, the resultant EOSs are referred 
to as the OI-EOSs. 

The OI-EOSs adopted here have a range of y < —200 
MeV fm 3 and 180 < A'o < 360 MeV, which results in the 
range of L as 0 < L < 180 MeV. Note that not only does 
such parameter range equally well reproduce empirical 
mass and radius data for stable nuclei, but also effectively 
covers even extreme cases [31]. We also remark that ac¬ 
cording to comprehensive re-analysis of recent data on 
the giant monopole resonance ener gies , Kq should be in 
the range of 250 < K 0 < 315 MeV [34|, while the gener¬ 
ally accepted value of K 0 is in the range of K 0 = 230 ± 40 
MeV [3a|. That is, systematic errors in experimentally 
determining Ko are still likely to be large. It is thus 
reasonable that the OI-EOSs used here have 247 sets of 
{y,Ko), i.e., the combination of 13 different values of y 
(y = -200, -220, -250, -300, -350, -400, -500, -600, 
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FIG. 2: (Color online) The parameter So is plotted as a func¬ 
tion of L for the 247 OI-EOSs. 


-800, -1000, -1200, -1400, and -1800 MeV fm 3 ) and 
19 different values of Kq (A'o = 180, 190, 200, • ■ ■ 360 
MeV). For these OI-EOSs, the corresponding values of 77 
are calculated, which are shown in Fig. |T] as a function 
of L. From this figure, one can observe that the depen¬ 
dence of 77 on L is much stronger than that on Kq. This 
is partly because uncertainties in L are relatively large 
compared with those of Kq and partly because the power 
of L in 77 is larger than that of Kq. 

Additionally, for comparison, we show Sq and L for 
the 247 OI-EOSs in Fig. [2] A strong correlation be¬ 
tween So and L was pointed out by [3l| , and is consistent 
with the values of So and L obtained via fitting to ex¬ 
perimental data on nuclear masses and radii on la level 
with the nuclear energy density functional for Skyrrne 
type interaction [36]]. As compared with this correla¬ 
tion, the correlation between 77 and L is equally strong. 
Via simultaneous observations of masses and radii of low- 
mass neutron stars, therefore, constraints on 77 and thus 
L would be available to some extent. Finally, in Fig. 
Q] we also show the fitting to the data of the 247 OI- 

EOSs, i.e., 77 = 6.59 ( i Mcv )~^ 3 MeV, which corresponds 
to Kq = 286.8 MeV. 

Constraints on 77 and L. Following the finding of 77 , 
we here give a possible constraint on 77 from observa¬ 
tions of masses and radii of low-mass neutron stars. Un¬ 
fortunately, however, no firm observational evidence for 
the presence of less-than-lMg neutron stars is available. 
Even more challenging is simultaneous mass and radius 
determination of low-mass neutron stars. Thermonuclear 
X-ray bursts in low-mass X-ray binaries help to deter¬ 
mine the masses and radii of the bursting neutron stars, 
although there exist many uncertainties both in theoreti¬ 
cal models and in observations. In fact, it is not straight¬ 
forward to determine the exact moment when the lumi¬ 
nosity reaches the Eddington limit at the star’s surface, 
if data for the photospheric radius expansion bursts are 
adopted to determine the star’s mass and radius. Addi¬ 
tionally, the color-correction factor defined as the ratio of 
the color temperature to the effective temperature of the 


source object is sensitive to the flux during the cooling 
tail as well as the model of neutron star atmospheres [29] . 
To minimize uncertainties in the theoretical models that 
determine the Eddington luminosities during the burst 
phenomena, Suleimanov et al. suggested using informa¬ 
tion from the whole cooling track in the X-ray bursts, 
and succeeded in obtaining the constraint on the mass 
and radius of the X-ray burster 4U 1724-307 located in 
the globular cluster Tarzan 2 by using various atmosphere 
models j4j. In this paper, we adopt their results, which 
imply a relatively low-mass neutron star, to obtain a con¬ 
straint on 77 and L. 

In particular, Suleimanov et al. adopted three atmo¬ 
sphere models with different chemical compositions, i.e., 
pure hydrogen, pure helium, and the solar H/He compo¬ 
sition with sub solar metal abundance Z = 0.3 Zq appro¬ 
priate for Terzan 2 [37]] • Then, assuming a flat distribu¬ 
tion of the distance from the Earth between 5.3 and 7.7 
kpc with Gaussian tails of lcr = 0.6 kpc, they obtained 
such constraints in the mass-radius relation within 90% 
confidence level as shown in Fig. [3J where the check¬ 
ered, filled, and shaded regions correspond to the con¬ 
stants obtained with the atmosphere models composed 
of pure hydrogen, pure helium, and the solar H/He with 
Z = 0.3-Zq, respectively. In addition to their results, we 
show the region ruled out by the causality, which is given 
by R < 2.824 GM/c 2 [H]. From this figure, one can ob¬ 
serve that the radius of the X-ray burster 4U 1724-307 
should be relatively large if the star’s mass has a canon¬ 
ical value of order 1.4 Mq. 

In Fig. [3] we also plot the stellar models constructed 
with several sets of the OI-EOSs and the Shen EOS [39j, 
where the corresponding value of 77 is written on each 
EOS. Here, we particularly focus on the stellar models 
for p c < 2.0po to avoid uncertainties in the EOS at high 
density due to the possible appearance of non-nucleonic 
components and/or the profoundness of three-neutron in¬ 
teractions as mentioned above. In Fig. [3l therefore, the 
upper end of each line corresponds to the stellar model 
constructed with p c = 2.0poJ for reference, we also show 
the stellar model for p c = 1.5po by putting a mark on 
each line. 

Now, assuming that the X-ray burster 4U 1724-307 has 
a canonical neutron star mass and that the EOS is univer¬ 
sal in the sense that all neutron stars can be constructed 
with a single EOS, one can conclud e from Fig.[3]that 77 is 
larger than ~ 130 MeV. Via L = \Jrf > /K q : 77 > 130 MeV 
leads to L > 110 MeV for Kq = 180 MeV, L >98 MeV 
for Kq = 230 MeV, L > 89 MeV for Kq = 280 MeV, and 
L > 78 MeV for Kq = 360 MeV. One can more clearly 
see the allowed region in the parameter space in Fig. [4l 
i.e., the region above the solid line, where we show the 
lines for 77 = 120 MeV (dashed line) and 140 MeV (dot¬ 
ted line) for reference. We remark that the density at 
the core-crust boundary strongly depends on the value 
of L II, which is a crucial property to determine the 
crust mass and moment of inertia. Combining Fig. 5 in 
Ref. [ 32 ] with the constraint on L obtained here from the 
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FIG. 3: (Color online) Allowed regions in the mass-radius re¬ 
lation obtained from the observation of the X-ray burster 4U 
1724-307 by Suleimanov et al. [lj, where they adopted three 
different atmosphere models, i.e., pure hydrogen (checkered 
region), pure helium (filled region), and the solar ratio of 
H/Hc with sub solar metal abundance Z = 0.3Zq (shaded 
region). On the other hand, the lines with marks denote the 
stellar models constructed from nine EOSs with different val¬ 
ues of tj (attached numbers) for p c < 2 . 0 po, where each mark 
corresponds to the mass and radius of a star with p c = 1.5po- 
Additionally, the upper left region is ruled out by the causal¬ 
ity [ 1 |. 


neutron star observations, the nucleon number density at 
the crust basis is expected to be around 0.07 fm -3 . 

Recently, a lower limit of observed neutron star radii 
has been additionally suggested from another object. 
That is, the neutron star radius in the low-mass X-ray 
binary 4U 1608-52 is predicted from the hard-state burst 
occurring during the low-accretion rate to be larger than 
13 km, if the neutron star mass is in the range of 1.2- 
2.4Mq j40|. This suggestion also indicates a large value 
of 77 , i.e., 17 > 100 MeV, which covers the constraint ob¬ 
tained from the X-ray burster 4U 1724-307. Thus, 17 is 
predicted to be larger than around 130 MeV from both 
of the astronomical observations. 

We conclude this section by noting that observations 
of neutron star oscillations could also tell us the prop¬ 
erties of neutron star matter. In fact, the gravitational 
waves emitted from oscillating neutron stars could pro¬ 
vide a possible way to see the neutron star properties, 
although they have not yet been observed directly. An¬ 
other possibility is the detection of electromagnetic waves 
associated with neutron star oscillations j4l], |42j . For¬ 
tunately, quasi-periodic oscillations have been discovered 
in the afterglow of giant flares observed from soft-gamma 
repeaters [16]. If such oscillations result from crustal tor¬ 
sional oscillations, a fairly strong constraint on L can be 
obtained by comparing the observed frequencies of quasi- 
periodic oscillations with theoretical predictions of the 
eigen-frequencies of the torsional modes |2(~il422f . In this 
way, we obtained two possible constraints on L. One is 
101.1 < L < 131.0 MeV, which explains all the observed 
frequencies lower than 100 Hz in terms of the crustal tor¬ 


FIG. 4: (Color online) Nuclear matter EOS parameters con¬ 
strained from 17 > 130 MeV, which correspond to the region 
above the solid line. For reference, 17 = 120 and 140 MeV 
are plotted with the dashed and dotted lines. We also dis¬ 
play the parameter space constrained from the observations 
of quasi-periodic oscillations in giant flares with the filled and 
checkered regions [ 22 ] (see text for details). 


sional oscillations, while the other is 58.0 < L < 85.3 
MeV if the second lowest frequency observed in SGR 
1806-20 would be excited by a different mechanism from 
the crustal oscillations. These constraints on L are also 
shown in Fig. [I] by the filled region for the former con¬ 
straint and by the checkered region for the latter one. 
As can be seen from this figure, the former constraint on 
L is more consistent with the constraint from 17 > 130 
MeV than the latter one. On the other hand, we remark 
that most of the terrestrial nuclear experiments suggest 
somewhat lower values of L 43|, although there still ex¬ 
ists large uncertainty in L [44j. We also remark that the 
EOS of pure neutron matter calculated within the chi¬ 
ral effective field theory favors smaller values of L (see 
[45l ] and references therein). Anyway, 17 (and L ) could be 
significantly smaller than our constraint, if constraints 
on M and R of several neutron stars (e.g., M) , other 
than the ones adopted in the present analysis, are taken 
for granted. 

Conclusion. We apply the classification of the EOS of 
neutron star matter in terms of 17 as developed in Ref. 
[ 27 } to the mass-radius relation constrained from observa¬ 
tions of the X-ray burster 4U 1724-307 and remark on a 
possible constraint on 77 , which gives important informa¬ 
tion on the density dependence of the symmetry energy. 
In order to obtain a better constraint on 77 , it would be 
significant to expect simultaneous mass and radius deter¬ 
mination of low-mass neutron stars from the neutron star 
interior composition explorer (NICER) by NASA and/or 
the large observatory for X-ray timing (LOFT) by ESA 
via observations of pulse profiles from hot rotating neu¬ 
tron stars [46j . 
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